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We propose an experimental realization of discrete quantum random walks using neutral atoms 
trapped in optical lattices. The random walk is taking place in position space and experimental 
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a quantum random walk to the classical random walk for increasing errors and decoherence. 
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I. INTRODUCTION 

The increasing effort to investigate theoretically and 
experimentally the possibility to construct and build an 
universal quantum computer is mainly motivated by the 
expectation that quantum computers offer a (possibly 
exponential) speed-up over classical computers. Despite 
the two celebrated milestones of Shor's factoring algo- 
rithm and Grover's database search algorithm — 
which both offer a speedup over their best (known) clas- 
sical counterpart — no constructive way to generate effi- 
cient quantum algorithms is currently known. One possi- 
ble direction of research is the adaption of known classical 
algorithms to the quantum mechanical case. 

Random walks on graphs play an essential role in var- 
ious fields of natural science, ranging from astronomy, 
solid-state physics, polymer chemistry and biology to 
mathematics and computer science g]. In particular, 
markov chain simulation has emerged as a powerful al- 
gorithmic tool and many classical algorithms are based 
on random walks. It is possible and hoped that quantum 
random walks allow in a similar way for a constructive 
search for new quantum algorithms. This justifies the 
increasing effort in the investigation of quantum random 
walks by several groups gl, || 0, |,|[ [l^, 0, H, • 
Different behavior of the quantum random walk — as 
compared to the classical one — have been reported un- 
der various circumstances. For instance, a very promis- 
ing feature of a quantum random walk on a hypercube, 
namely an exponentially faster hitting time as compared 
to a classical random walk, has been very recently found 
(numerically) by Yamasaki et al. and (analytically) 
by Kempe @ . 

In this paper, we consider the simplest and best- 
studied version of a random walk, namely the discrete 
Hadamard walk on a line or a circle, first studied by 
We propose an experimental implementation of the 
quantum random walk using neutral atoms trapped in an 
optical lattice. In contrast to the recently proposed im- 
plementations using ion traps put forward by Travaglione 
and Milburn [Il2| a nd microwave cavities put forward by 
Sanders et al. ||14| , in our proposal the random walk is 
taking place in position space and several hundred steps 
may be implementable even with present day technology. 



The paper is organized as follows. In Sec. II, we com- 
pare features of the classical and quantum random walk 



on the line and introduce some basic notation. Sec. Ill 
provides a description of the physical set-up using opti- 
cal lattices and implementations of the Hadamard walk 
on a line using this set-up. Possible imperfections and 
their influence on the quantum features of the walk are 
discussed in Sec. IV. We summarize and conclude in Sec. 
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II. 



CLASSICAL VS. QUANTUM RANDOM 
WALKS 

A. Classical random walk on a line 



Consider an infinite line with allowed (integer) posi- 
tions Xk = k, k ^ "L and a particle which is initially 
located at position a;o = 0. We consider a step- wise 
evolution in such a way that at each step, the par- 
ticle moves with probability 1/2 one step to the left, 
x{n) = x{n — 1) — 1, and with probability 1/2 to the 
right, x{n) — x{n — 1) + 1. After n steps, the probability 
Pciassicai(f^, k) to find the particle at position x^ is given 
by 
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Note that if n is even [odd] , only even [odd] positions are 
occupied. The standard deviation of the distribution is 
■\/n, which implies a spreading time proportional to ^/n- 
The probability to observe a particle at distance of order 
n from the origin decreases exponentially with n and is 
zero Vno, hq > n. 



B. Quantum random walk on a line 

A quantum mechanical analogy for the classical ran- 
dom walk would be a particle whose state evolves at each 
step into a coherent superposition of moving one step to 
the right and one step to the left. One readily finds Q 
that unitarity of the evolution implies that one has to 
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consider a particle with internal degrees of freedom to 
achieve this aim. We consider a particle with two inter- 
nal degrees of freedom, which can move on an infinite line 
with integer positions j^, The corresponding Hilbert 
space H = Hi ® Tix is given by Hi —(T^, the internal 
state of the particle, and Hx ~<I°°, the position space 
with basis states \k)x G Hx, —oo < k < oo, correspond- 
ing to the particle located at the fc**^ lattice site.. 

The internal state of the particle, |0),|1) determines 
the direction of the particle movement. If the internal 
state is |0), the particle moves to the left, while it moves 
to the right if the internal state is This operation is 
described by the unitary controUed-shift operation, 

S = |0)(0| ®Y.\^- l)(fc| + |1)(1| ® 5] |fc + l)(fc|, (2) 

fe k 

i.e. S\0)i (gi \k)x = \0)i ® \k ~ l)x and S\l)i ® \k)x ^ 
\l)i®\k + l)x. 
We also introduce the Hadamard operation, 

which acts on Hilbert space Hi such that H\0) = 
1/^2(10) + |1)), H\l) = l/\/2(|0) - |1)). The particle 
is initially prepared in state \ipo) — l/\/2(|0)/ + ® 
|0)x- Each step of the quantum random walk — which 
is also called Hadamard walk — consists of applying the 
Hadamard operation, H (gl, followed by the controUed- 
shift operation S. Let j-^n) ~ {SH)"\ipo) be the state 
of the system after n steps. The probability p{n, k) to 
observe a particle at position k after n steps is given by 

=trx(|A;)x(fc|tr/(|V'„)(^n|)), (4) 

and may be compared to the probability distribution 
Pciassicai(?^) of the classical random walk. The proba- 
bility distribution p(n) has been analyzed in detail in 
Ref. §. While 

Pclassical 

(n) is given by a binomial dis- 
tribution — which is for large n well approximated by 
a gaussian — , no such simple form exists for p{n). The 
standard deviation of the distribution pciassicai(fi) is \/n, 
while p{n) is almost uniformly distributed in the interval 
(— n/\/2, n/y/2) and the standard deviation is linear in n. 
This implies that the spreading time for a particle goes 
like ^/n for the classical random walk, while in the quan- 
tum random walk it scales linearly with n. This provides 
an essential different behavior of the quantum random 
walk that follows from the possibility of interference in 
the quantum mechanical case. 

In a similar way, the (quantum) random walk on a 
circle is defined using a position space Hx = d/^ with 
periodic boundary conditions, i.e. \k)x = \k'niod{N)) x 
for some finite N. Also in this case, a quadratic speed-up 
of the quantum random walk compared to the classical 
random walk is found in the spreading time of the particle 

. Walks on general graphs can be defined in a similar 
way [§. 



III. IMPLEMENTATION IN OPTICAL 
LATTICES 

In this section, we discuss possible implementations of 
the quantum random (Hadamard) walk on a line or on 
a circle using neutral atoms trapped in periodic optical 
potentials (for a review see e.g. [l7|l). 

A. Physical set— up 

We consider two identical one-dimensional optical lat- 
tices, each of them trapping one of the internal states 
|0), |1) of a neutral atom. For example, one may use alkali 
atoms with a nuclear spin equal to 3/2 (*^Rb, ^'^Na), and 
choose the hyperfine structure states |F = l,m/ = 1) 
[\F = 2, m/ = 2)] to represent |0)[|1)] respectively. 

Each lattice consists of a periodic optical potential 
with period d. The optical potentials are formed by the 
standing waves resulting from two counter-propagating 
traveling waves with the electric fields forming an an- 
gle of 29, the so called linZlin configuration. By chang- 
ing 9, the right and left circular polarized components 
CT^ of the standing waves forming the total electric field 
can be shifted with respect to each other, E'^{x,t) — 
i?oe^*'^*[e+ sin(fca; + 9) + sin(fca; — 9)]. We have de- 
noted k = v/c the laser wavevector, Eq the amplitude 
and e± unit right and left circular polarized vectors. The 
potentials "seen" by the atoms in the internal states 
|0),|1) s.reVo{x,9) = [Kn,=i/2(x, 0) + 3K™^=_i/2(x, 0)]/4 
and Vi{x,9) = ^^,=1/2(2;, 6*), where ^^^^±1/2(2;, 6*) = 
a\Eo\'^si^^{kx±9) @. 

This basic architecture can be used for quantum state 
control of neutral atoms in optical lattices and it con- 
stitutes the basis of the proposals for quantum computa- 
tion in such systems |l^ . As in the proposal of Ref. 
jlSf , we make use of the fact that a relative movement 
of the two lattices, i.e. the trapping potentials Vo,Vi, 
can be achieved by varying the angle 9. In particular, 
starting with = 0, the respective minima of the poten- 
tials Vb, Vi coincide and by changing 9 from to 7r/2, the 
potentials Vo,Vi move in opposite directions until their 
respective minima coincide again. Note that the shape 
of the potential Vb changes as it moves. 

B. Implementation of the Hadamard walk 

We consider a single neutral atom at position ccq = 
and the case where lattice — which traps the internal 
state |0) of the neutral atom — moves with constant ve- 
locity to the left, vq ~ — w, while lattice 1 — which traps 
the internal state |1) of the atom — moves with constant 
velocity vi — v io the right. The initial position of the 
lattices is such that the minimum of a potential well is 
located at position ccq at to = 0- The lattice movements 
are used to implement the controUed-shift operation (see 
Eq. (ph), while laser pulses allow one to manipulate the 
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internal state of the atom and thus to select the the corre- 
sponding trapping potential (and therefore the direction 
of the movement). 

Given that the atom is initially prepared in state 
l/-\/2(|0) + at position xq — 0, the application of 
the Hadamard operation (see Eq. ^) to the internal 
state of the atom at times t„ = nd/v readily implements 
the quantum random walk on a line using this set-up. 
The spatial probability distribution of the atom at time 
tn — i.e. the probability to observe an atom at posi- 
tion kd, —n < k < n sX time tn — corresponds exactly 
to p{n,k) (Eq. (^)) of the one-dimensional Hadamard 
walk after n steps. A simple fluorescence measurement 
— together with several repetitions of the experiment — 
allows one to measure this distribution. 

To justify this statement, note that a single atom ini- 
tially at position ccq = which is prepared in state |0) 
moves with constant velocity to the left, x{t) — 
Xq — vt [right, x(t) = Xq + vt] respectively. After a time 
tn = nd/v, the position of an atom is shifted by exactly 
n lattice periods nd and the two lattices are again on 
top of each other. By changing the internal state of the 
atom e.g. at time tn from |0) to one can switch 
between the corresponding trapping potential and thus 
change the direction of the movement. Note that it is 
important to make such changes of the internal state of 
the atom only when the two lattices are on top of each 
other, to ensure that the atom remains trapped in one 
of the potentials. Coherent superpositions of two inter- 
nal states behave likewise. In the case of Rubidium with 
|0) = \F = 1,TO/ = 1), |1) ^ \F = 2,m/ = 2), one can 
use standard Raman pulse or microwave techniques to 
realize the Hadamard rotation by using fast laser pulses. 

Note that on the Blochs sphere, the Hadamard opera- 
tion corresponds to a rotation of an angle tt around the 
axis u — l-\/2(ei + e^). This corresponds to a 7r~pulse 
rather than a 7r/2-pulse in the usual terminology of quan- 
tum optics. Up to an irrelevant global phase, one may 
also achieve the Hadamard operation by a sequence of 
three 7r/2-pulses, H cx e-''^^'^''- e'^'^^^'^y e"'"'/^'^- . Exper- 
imentally it may be easier to use a 7r/2-pulse correspond- 
ing to the transformation — e^'^/'*'^=" instead of the 
Hadamard operation and prepare the atom initially in 
state l/V2{\0) + |1)). This also leads to a symmetric 
probability distribution for all times t equivalent to the 
one resulting from the standard Hadamard walk. 

The Hadamard operation has to be applied at all lat- 
tice sites, which can be easily achieved by using a non- 
focused laser beam. In fact, such a homogeneous oper- 
ation is much easier to implement than individual 
operations on specific lattice sites. This is due to the fact 
that in current experiments, the lattice period d ~ 425nm 
— which is limited by the optical wavelength — is smaller 
than the best achievable focusing width of the laser 
beams, w ~ Ifim p3| . In the fluorescence measurement, 
one can either detect unselectively both internal states 
|0) and |1) to reveal information about the position of 
the atom, or one may use selective resonance fluorescence 



methods. In the latter case, addition application of a ran- 
dom ax operation (7r-pulse) before the measurement is 
required to remove the dependence on the internal state 
of the spatial probability distribution. Provided the atom 
was initially prepared in state l/\/2(|0) -f the prob- 
ability distribution is symmetric when tracing out the 
internal state of the atom, however the distributions con- 
ditioned on the internal state of the atom are asymmetric 
and have mirror symmetry, which explains the additional 
application of a random operation. Notice the phase 
i in the initial state, which is important to ensure sym- 
metric behavior of the random walk. 

We would like to emphasize that the procedure 
sketched above is readily implementable with existing 
technology. It does not require addressability of indi- 
vidual lattice sites. 

The essential requirements are that the internal states 
of the atom — as well as their coherences — are sufficiently 
stable and that the particle remains trapped in the po- 
tential throughout the procedure. This can be satisfied 
if the movement of the lattice is sufficiently slow [that 
is f <C i^osc, where Vosc ~ iquj is the rms velocity of 
the atoms in the vibrational ground state, oj is the ex- 
citation frequency and on is the size of the ground state 
of the trap potential [Q] such that the atom stays in 
the ground state of the potential during motion. This 
condition can be relaxed, as will be discussed in the sub- 
sequent section, and one can also allow non-adiabatic 
velocity profiles. The coherence of the internal state is 
mainly affected by fluctuations the intensity and phase of 
the trapping lasers |2j] as well as magnetic field fluctua- 
tions which may lead to uncontrollable energy splittings 
between the internal levels. We will address some of these 
issues in the next section. Given that these noise sources 
can be controlled sufficiently well, the number of steps of 
the random walk one can perform is only limited by the 
spontaneous emission lifetime of the atom in the lattice, 
which is at the order of several seconds. This corresponds 
to a maximum number of about n ~ 10* time steps, as- 
suming t\ sa 100/iS— 1ms which respects the adiabaticity 
requirement for lattice shifts. Note that the implementa- 
tion of several hundred time steps of the Hadamard walk 
corresponds to a spatial width of the quantum distribu- 
tion at the order of millimeters. 



Improved implementation of the Hadamard 
walk 



From a practical point of view, there are a number 
of difficulties with the procedure proposed in Sec. Ill B . 
For example, the laser pulses to implement the Hadamard 
rotation have to be fast compared to the timescale of the 
lattice movement. In addition, if the internal state of an 
atom is changed e.g. from |0) to |1), this imphes a sudden 
momentum change of the atom, as it is no longer trapped 
in the left-moving lattice but in the right-moving one. 
This momentum change may lead to heating of the atom. 
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and the atom may eventually even escape from the trap. 

Another practical difficulty one faces in current exper- 
iments is concerned with dephasing of the internal states 
of the atom. In particular, uncontrollable time and space 
dependent magnetic fields lead to energetic shifts of the 
internal levels, which result in relative phase shifts de- 
stroying the coherence of the system . 

In this section, we propose a slight modification of the 
implementation suggested in Sec. IIIB. This scheme 



is based on symmetrizing the procedure and avoids the 
problems mentioned above. Instead of moving the lat- 
tices with constant velocity, they oscillate around the 
central position xq — 0. In the simplest case, the move- 
ment of the lattices is harmonic and may be described as 
follows 



x'{t) = d/2{cosujLt- 1), 
x"{t) = -d/2{cosL0Lt - 1), 



(5) 



where x'{t) [x"{t)] is the position of the central poten- 
tial well of lattice [1] respectively. The oscillation fre- 
quency lol is chosen in such a way that the adiabatic- 
ity requirement for the lattice movement — i.e. that the 
atom remains in the motional ground state throughout 
the procedure — are well fulfilled, which leads to oscilla- 
tion times at the order of 100/iS to ms. Note that one 
may replace the simple harmonic movement of the lat- 
tices by more complicated profiles, either specially suited 
to meet adiabaticity requirements such as the one pro- 
posed in Ref. jl^, or specially designed in such a way 
that adiabaticity requirements need not be matched, but 
are replaced by the weaker condition that the atoms — 
after moving the lattice b y o ne period — are again in 
the motional ground state |Q . These specially designed 
profiles may allow for movement times at the order of 
a few tens of microseconds for a shift of one lattice pe- 
riod. The profiles of Ref. need to be adopted 
in such a way that the lattices oscillate around a central 
position, which may be accomplished by choosing the 
original profile until the lattice is displaced by one lattice 
period and the velocity is zero to good approximation, 
and then use the time-inverse profile. In this way, pe- 
riodic lattice movements are readily achieved. In what 
follows, we will restrict our discussion to harmonic lat- 
tice movements, which is sufficient to illustrate the ideas 
of the improved implementation of the Hadamard walk. 
This is realized if, in addition, ax operations (7r-pulses) 
are applied at times t„ = mr/ujL,n = 1,2,3,... to all 
lattice sites, with the effect |0) |1). 

Under these conditions, an atom initially in state |0) is 
trapped in lattice and starts moving to the left together 
with lattice 0, where at time ti it is located at position 
~d and the internal state changes to |1) due to the ap- 
plication of ax- The particle is therefore now trapped in 
lattice 1, where which moves to the left in the interval 
(ti, ^2), leaving the trapped atom at position —2c? at time 
t2 etc.. On the one hand, we have that at times tn — 
when manipulations of the internal states of the atoms 
are performed — the two lattices are on top of each other 



and the velocity of the atoms is zero, which overcomes the 
first difficulty mentioned above. One could in principle 
also stop the lattice movement at these times until ma- 
nipulation of the internal states is achieved, which allows 
to drop the requirement that manipulation of the inter- 
nal states of the atom have to be fast compared to the 
timescale of the lattice movement. On the other hand, 
since within a time span of 2ti the atom is both for time 
ti in state |0) and for ti in state |1), relative phase shifts 
between internal states |0) and |1) become an irrelevant 
global phase shift. Also fluctuation of magnetic fields 
become irrelevant, provided the timescale of the fluctu- 
ations is much larger than 2ti and the spatial variation 
is negligible within 2c? ]2^ . These requirements are well 
fulfilled e.g. for the typical 50IIz background noise and 
oscillating times ti w 100/is— 1ms when using harmonic 
lattice movements. 

The implementation of the Hadamard walk using this 
set-up is straightforward. After the ax operations at 
times tn, the Hadamard operation H (see Eq. (||)) is ap- 
plied at times tn if n even, while H' = axHa^ is applied 
at times tn if n odd. That is, at to H is applied, while 
at tn H'ax = (JxH [Hax\ is applied if n is odd [even]. 
The use of the operation H' instead of H results from 
the interchanged role of the internal states |0) and |1) for 
even/odd n. One readily checks that in this way, after 
time t„, n steps of the quantum random Hadamard walk 
are implemented, provided the internal state of the atom 
is l/v^(|0) + i|l)). The additional application of H,H' 
does not infiuence the symmetrization discussed in the 
previous paragraph. Fluoresce nce m easurements can be 
performed as described in Sec. IIIB. 



D. Bounded random walks and random walk on a 
circle 

Using above set-up one can also implement an one- or 
two-side bounded random walk on a line [|, |l|. Such a 
bounded random walk is similar to its unbounded coun- 
terpart, however at a certain locations xi, xi barriers are 
introduced and the walks ends once a particle reaches one 
of these locations. A one side bounded random walk may 
e.g. contain a barrier at X\ = 100, while xi —00. Such 
a barrier can in the optical lattice set-up for instance be 
implemented by shining a laser at a certain location xi, 
which couples both state |0) and |1) to a fast decaying 
auxiliary level. In such a way, a (position) measurement 
of the atom projecting onto Pi = is performed. 

A modification of the trapping geometry may also al- 
low for the implementation of a random walk on a circle, 
following the ideas of a recent proposal by Burke et al. 
|^l| , using an evanescent field of a linear waveguide and 
a ring resonator for trapping and guiding atoms (see also 
Ref. [Q). In this way, the periodic trapping potential 
can be modified such that lattices sites are located on a 
circle, forming a regularly spaced pattern. Movements 
of the trapping potentials result in this case in a circu- 
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lar movement of the lattice sites and thus of the trapped 
particle. By means of Hadamard rotations together with 
lattice movements, a random walk on a circle could be 
implemented. Measurements and manipulations of the 
trapped atom can b e perfor med in the same way as dis- 
cussed in Sec. Ill B and Sec. QIC. Note that in contrast 
to the quantum random walk on the line, the random 
walk on the circle can not be implemented using exist- 
ing experimental set-ups but rather relies on a proposed 
scheme. 



E. Using 2D set— up to measure probability 
distribution of a ID random walk 

In current experiments, two or three dimensional lat- 
tice arrays are used rather than ID-arrays. In this 
case, four [six] interfering laser beams constitute the two 
[three] dimensional trapping potential. One can make use 
of such a set-up to directly measure the probability dis- 
tribution of the one-dimensional random walk. Consider 
a two dimensional lattice, which may be loaded from a 
BEC [Q. We assume that the 2D lattice is loaded in 
such a way that in one dimension, say x, only the central 
lattice site xq is occupied, while in the other dimension, 
say y, all lattice sites are occupied in a regular way with 
one atom per site. This can be accomplished, for exam- 
ple, by first realizing a Mott transition in a 3D optical 
lattice from a BEC to a Mott insulator state as in re- 
cent experiments of the Munich group |2^, following a 
theoretical proposal of Ref. [^. In a second step, one 
can e.g. use of a gradient magnetic field to selectively 
address specific lattice layers and to deplete unneeded 
lattice sites 23 . One ma y a lso apply methods similar 
to the ones used in Ref. |28|| to achieve uniform filling 
factors. 

Such a configuration allows for a parallel sampling of 
one dimensional random walks (in the x direction), by 
moving the lattices in the x direction only and apply- 
ing H adamard r otatio ns to all atoms as described in Sec. 



IIIB and Sec. [II C, Each of the atoms at position 
Uk = kd independently performs a random walk. A flu- 
orescence measurement, e.g. a projective picture along 
the y-axis, allows to directly measure the corresponding 
probability distribution p{n) provided the number of lat- 
tice sites in y direction is big enough. Otherwise, the 
required number of repetitions of the experiment to de- 
termine the probability distribution is decreased. 



IV. EXPERIMENTAL IMPERFECTIONS 

Although existing experiments using optical lattice 
systems ofi^er high accuracy in both coherent storage and 
manipulation of the atoms, different kinds of errors may 
influence the ideal evolution. These errors may disturb 
or even destroy typical features of the quantum random 
walk, such as linear spreading time. In this section, we 



concentrate mainly on errors in the coherence of the in- 
ternal states of the atoms. We thereby observe a tran- 
sition from quantum mechanical to classical behavior of 
the random walk for increasing errors. We expect these 
kind of errors to constitute the dominating part in the 
experimental imperfections. On the one hand, errors in 
lattice movements may lead to motional excitations of the 
atom. If sufficiently small these should however not effect 
the essential behavior of the system. On the other hand, 
the internal states of the atom are influenced by decoher- 
ence resulting from e.g. uncontrollable phase shifts, im- 
perfections in the manipulation by means of laser pulses 
as well as fluctuations in the trapping potential during 
lattice shifts. One may distinguish between errors intro- 
duced by manipulations of internal state of the atom and 
errors in the position space of the atom, e.g. introduced 
by tunneling of atoms between neighboring lattice sites. 
While the former always keeps the structure of the ideal 
Hadamard walk that after n time steps, only even [odd] 
lattice sites are occupied if n is even [odd] respectively, 
errors in the latter lead to occupations of all lattice sites. 

We use two simple models to investigate imperfections 
in the coherence of the internal state of the atom. We 
treat decoherence effects and errors in the manipulation 
of the atom (imperfect Hadamard operations), i.e. op- 
erations acting on the total Hilbert space Hi (Xi Hx as 
[// (8" Ijjf, in a joint way. In the first model, we assume 
that the desired manipulation of the internal states of 
the atom, U, is performed with probability p at each 
time step, while with probability 1 — p a completely de- 
polarized, random state is produced. The parameter p 
serves not only as a measure of the accuracy of the oper- 
ation — where p = 1 describes perfect operations, while 
p = corresponds to a completely random operation — 
but also includes other decoherence eS'ects due to storage 
errors or phase fluctuations as well as lattice movements. 
Such a covariant error model reflects our limited knowl- 
edge about the specific type of error which occurred in 
the system. This error model has also been used in other 
contexts and is described by the following mapping 



£{p) - pUipUj + (1 - p)l/2h (g> tr/(p). 



(6) 



Note that this model is equivalent to a (par- 
tially) depolarizing channel, £{p) = pUipUj -f (1 — 

P)^/^Y^'k=Q '^li "' P'^'k ^ ^ where ak are Pauli matrices with 
ao = 1. 

The second model only includes phase errors and is mo- 
tivated by the expectation that phase fiuctuation may be 
the dominating part of errors occurring in optical lattice 
systems. This error model is described by the following 
mapping 



£'ip)^p'UipUj + il-p')Uia., 



3 PCTg 



(7) 



Note that if the optical potentials are not very deep, 
tunneling between neighboring sites may occur as well. 
We have used a simple model of incoherent tunneling - 
affecting only the position of the atom — which is given 
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by the following mapping |Q 

£"{p) = + (1 - q)/2{U+pUl + U-pUi), 



(8) 



where U± = Ij ^ J2k 1^ ^ the unitary shift 

operator which moves the particle either one position to 
the left, U- or right, U+. That is, with probability q 
nothing happens — and thus in total the desired evolution 
occurs — , while with probability (l — q) a tunneling of the 
atom to one of the neighboring lattice sites occurs. 

We have performed numerical simulations to investi- 
gate the influence of these kinds of errors on the quan- 
tum random walk on the line, where we first assumed 
that errors affect only the internal state of the atom. 
Figure |l] is based on error model 1 and shows the prob- 
ability distribution after n = 200 steps of the random 
walk for different error parameters. Note that for com- 
pletely random operations, i.e. parameter p = 0, the par- 
ticle performs exactly a classical random walk. This can 
easily be understood by observing that an internal state 
1/21 = 1/2(|0)(0| + |1)(1|), when applying the controUed- 
shift operation (i.e. the lattice movement), has the effect 
that the particle moves with probability 1/2 to the left, 
while with probability 1/2 it moves to the right. In con- 
trast to the quantum random walk, the resulting state 
is an incoherent superposition of the two possible states, 
which can be described classically and thus no interfer- 
ence effects (as in the quantum random walk) occur. The 
internal state plays the role of a classical coin. One ob- 
serves from Fig. that with increasing errors (decreasing 
parameter p), the probability distribution changes from 
the quantum mechanical one to the classical one. Even 
for errors at the order of several percent, typical quan- 
tum mechanical features of the probability distribution 
after a few hundred steps are clearly visible, in particular 
occupation in the interval {y/n,n/ \/2) can be observed. 
A similar simulation was performed using error model 
2 (Eq. 1^). The observed behavior of the system under 
this kind of error is very similar to the one shown in 
Fig. 1^ using error model 1. Fig || shows the probability 
distribution after different number of steps of the ideal 
[imperfect] quantum random walk assuming only phase 
errors. 

In the following wc assume both internal and external 
errors (tunneling), described by Eq. ^ and Eq. (||) 
respectively. As shown in Fig. ||, the essential effect of 
incoherent tunneling is that the probability distribution 
is smeared out. 

We have also considered a one-side bounded random 
walk. Fig. ^ shows the probability to observe the atom 



at the barrier at position a; = —10 as a function of the 
number of steps. As expected, the exit probability for 
the quantum random walk is smaller than 1, while it 
approaches unity for the classical walk ||l^ . 

Although no reliable estimates for the parameters p in- 
cluding all possible imperfections and decoherence effects 
are available, errors at the order of several percent are 
still tolerable to observe a clear quantum behavior of the 
random walk even after a few hundred steps. This seems 
to be experimentally achievable. In turn, the distribution 
measured in the experiment can be used to determine the 
degree of coherence of the system, in particular the qual- 
ity of the implemented operations. This may also serve 
as a test on the suitability of optical lattice systems to 
perform general purpose quantum computation, follow- 
ing the proposals of Ref. ||lq| and Ref. |3l|. 



V. SUMMARY AND CONCLUSIONS 

We have proposed to use neutral atoms trapped in op- 
tical lattices to implement quantum random walks on 
the line and on the circle. The random walk is per- 
formed in position space by periodically shifting the lat- 
tices and manipulating the internal states of the atom(s) 
by homogeneous laser pulses. Read-out of the result- 
ing probability distribution is performed via fluorescence 
measurements. Due to long life-times of the trapped 
atoms and efficient manipulation techniques, experimen- 
tal realizability is expected with present-day technology. 
We have also investigated the influence of decoherence 
and imperfections in manipulation of internal state of 
the atoms and showed a transition taking place from the 
ideal quantum random walk to the classical random walk 
for increasing errors. Errors at the order of percent seem 
tolerable to still observe a clear quantum behavior of the 
walk after a few hundred steps. 
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position are occupied and plotted, which explains the lower 
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